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THE AVERAGE LENGTH OF A TRAJECTORY IN A CERTAIN 
"rn" 1 BILLIARD IN A FLAT TWO-TORUS 

o 
o 

C\J ' F.P. BOCA, R.N. GOLOGAN AND A. ZAHARESCU 

o 

<D ■ 

Q Abstract. We remove a small disc of radius e > from the flat torus T and consider a 
point-like particle that starts moving from the center of the disk with linear trajectory under 
' angle uj. Let t e (uj) denote the first exit time of the particle. For any interval I C [0,27r), 

any r > 0, and any 8 > 0, we estimate the moments of f e on / and prove the asymptotic 
formula 

Jr r E {uj)duj = c r \I\E- r + Os(e' r+ i- S ) as e -> 0+, 

■ where c r is the constant 

-S : 12 T/ , _! i - (i - .t)^ 1 \ 

■ J \ rx(l - x) (r + — x) ) 

o 

A similar estimate is obtained for the moments of the number of reflections in the side 
cushions when T 2 is identified with [0, l) 2 . 

> 
oo 
O 

CSJ ■ 1. Introduction and main results 

O i 

For each < e < ^ we consider the region 

Z £ = {z £ M 2 ; dist(z,Z 2 ) > e} 

^ ' and the first exit time (also called free path length by some authors) 

T e (z, oj) = inf{T > ; z + tuj £ dZ £ }, z £ Z E , uj £ T, 

of a point-like particle which starts moving from the point z with linear trajectory, velocity 
uj, and constant speed equal to 1. This is the model of the periodic two-dimensional Lorentz 
gas, intensively studied during the last decades (see g], |7j, [S], 0, [TO], [IE], [12], U2], d, 
CS],!IE],!II],[IH],!2n],!2I],!2i,!3I],!S2]fora non-exhaustive list of references). The phase 
space of the system consists in the range of the initial position and velocity and is one of the 
spaces Y e xT with the normalized Lebesgue measure, or T,f = {(x, y) G 8Y £ x T ; uj ■ n x > 0} 
with the normalized Liouville measure. 

Equivalently, one can consider the billiard table Y £ = Z £ /1? obtained by removing pockets 
of the form of quarters of a circle of radius e from the corners. The reflections in the side 
cushions are specular and the motion ends when the point-like particle reaches one of the 
pockets at the corners. In this setting t £ (z,lo) coincides with the exit time from the table 
(see Figure P). 
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This paper considers the situation where the trajectory starts at the origin O = (0,0). In 
this case the phase space only consists in the range of the initial velocity of the particle. It 
is given by the one-dimensional torus T and can be reduced, for obvious symmetry reasons, 
to the interval [0, f ] . From the point of view of Diophantine approximation this corresponds 
to a homogeneous problem. We shall be concerned with estimating the moments of the first 
exit time f e (uj) = t £ (0,uj) as e — > + when the phase space is the range [0, f] of the velocity 
uj. This question was raised by Ya. G. Sinai in a seminar at the Moskow University in 1981. 
We answer the question by supplying asymptotic formulas with explicit main term and error 
for all the moments of f F in short intervals as follows: 



Theorem 1.1. For any interval I C [0, ?] and any r, 5 > 0, one has 

'0 r) s(el- s ) if r ^2 
O r $(e* - *) ifr = 2 



e r j T r e {u)dw = c r \I\ + { ^ ' ) i_ s J _ c ^ n ase^0 + , 



where 



1/2 

nil. r-i - i-d-xY 1 - ( I 







x(x 



■K 2 J \ rx(l — x) (r + l)x(l — x 



The mean free path length is in this case 

l 

f e (oj) duo 



tt/4 

4 f _ . , , ci 12 In 2 0.421383 

- / T £ [U)) duj r = . ~ , 

TT £ TT Z 2e £ 



1/2 

lim c r = ft / — t r dx 



Note also that 



r-»0+ ' IT 2 J x(l 
o 

To prove Theorem 11.11 we first replace the circular scatterers by cross-like scatterers [m — 
£, m + e] x {n} U {m} x [n — e, n + e], m, n G I? \ {(0, 0)}. 1 We denote by l e (uj) the free path 
length in this situation, and first prove 

Theorem 1.2. For any interval I C [0, j] and any r,a,S > 0, one has 

-2a 

t * 

Os{£ 



r [jr, w f dx , \Or,5{^~ 2a ~ 5 + W a ) $ T ± 2 

e / l'(uj)duj = c r / h< i j as e->0 + . 

7 eV ; 7 cos r x \o s (£2- s ) ifr = 2 

We consider the probability measures p^. and [i\ on [0,oo), defined by 

filtf) = 4r / /KH) ^, = 4r y /KH) dw, / e C c ([0, oo)). 



1 Actually it is not hard to see that for u) £ [0, j] the result for cross-like scatterers is asymptotically the 
same as when using vertical scatterers {m} x [n — e, n + e]. 
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Their supports are all contained in [0, y/2] as a result of Lemma ETT1 Moreover, we infer from 
Theorems II. II and II. 21 that their moments of order n £ N* are of the form 2 

= -j f e » du = c n + j } - O n ,s(es- 5 ); 
I 

I I 

These asymptotic formulas show in particular that jll(X n ) and yu^(X") converge to the main 
terms as e — > + . The Banach-Alaoglu and Stone- Weierstrass theorems now lead to 

Corollary 1.3. There exist probability measures fl and fi 1 on [0, y/2] such that 

fil — > p, and fj, e —>■ fi 1 weakly as e — > + . 

Moreover, the moments of ft, and [i 1 are 

oo 

J t n djx{t) = c n 


and respectively 

oo 




Figure 1. The trajectory of the billiard 

Besides, we estimate the average of the number of reflections R e (uj) in the side cushions of 
the billiard table in the case of circular scatterers and prove 

Theorem 1.4. For any interval I C [0, j] and any r, 5 > 0, one has 

e r j R r e {uj) duj = c r j {smx + cos x) r dx + O r ^{e&~ s ) as e — > + . 
I I 

2 We denote N = {0, 1, 2, . . . } and N* = {1, 2, 3, . . . }. 



4 



BOCA, GOLOGAN AND ZAHARESCU 



Again, we first consider the case of cross-like (or vertical) scatterers, let R e {u) denote the 
number of reflections in the side cushions of the billiard table in this case, and prove 

Theorem 1.5. For any interval I C [0, j] and any r,a,5 > 0, one has 

e r J R r £ (co) alto = c r J (l + t&nx) r dx + O rtS (ei- 2a - s + \I\e a ) ase^0 + . 
I I 

We may also consider the probability measures u £ and v l e on [0, oo) associated with the 
random variables eR e and eR £ , and defined by 

*/(/) = 4r / f{eRs{u)) du, vlU) = f[sReH) dw, f G C c ([0,oo)). 

I I 
From Theorems 11,41 and 1 1 1 51 we derive 

Corollary 1.6. There exist probability measures v 1 and v l on [0, y/2] such that 

9 £ — > v 1 and 1,1 "' 
Moreover, the moments of v 1 and v 1 are 



v £ — ► v and u e —> v as e —> + . 



and respectively 



J t n dD\t) = %J (sinx + cosx) n dx, 
o I 

oo 

J t n dv\t) = jj- J (1 + tan x) n dx. 



In the case I C [j,^] one gets formulas similar to the ones in Theorems ll.il H~2l 11.41 and 
ll.5l after performing a symmetry with respect to a diagonal of the square, i.e. replacing (a, (3) 
by (§-/?,§ -a). 

The proofs make use of techniques employed in the study of the spacings between Farey 
fractions, pioneered in |2H> 123- and furthered recently in [3J, [1], pQ, [2E] where estimates 
for Kloosterman sums are being used. The first step consists in proving Theorems 11.21 and 
11.51 which refer to the case of cross-like or vertical scatterers. In this case one can directly 
take advantage of the fact that the intervals I a j q = [^p, s ^], with | Farey fraction of order 

Q = [-], provide a covering of [0, 1] such that two intervals I a / q and I a i i q i overlap if and only 

if | and ^ are consecutive Farey fractions of order Q. 

Finally, the case of circular scatterers is settled by partitioning the range I into 
intervals of equal size for a convenient value of the exponent 6, and replacing the small circles 
of radius e first by vertical scatterers of type {m} x [n — e_(m, n),n + £+(m, n)], and finally 
by scatterers of type {m} x [n — e, n + e] for appropriate choices of e±(m, n) and e. 

It should be possible in theory to compute the densities of the limit measures from their 
moments using either the Cauchy transform or the inverse Mellin transform. An attempt of 
this kind does not seem to easily lead however to a tractable formula for these densities. The 
convergence of the measures \Ju\ and v\ was proved in a different way and the limit measures 
were explicitly computed in 0. 
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Techniques using Farey fractions and Kloosterman sums were recently used in [0] to es- 
tablish the existence, and compute the distribution, of the free path length for the periodic 
two-dimensional Lorentz gas in the small-scatterer limit in the case where the trajectory does 
not necessary start from the origin, and one averages over both initial position and initial 
velocity. 



This is the final version of the paper with the same title, circulated as preprint math. NT/01 10208 



2. Farey fractions and Kloosterman sums 



For each integer Q > 1, let Tq denote the set of Farey fractions of order Q, i.e. irreducible 
fractions in the interval (0, 1] with denominator < Q. The number of Farey fractions of order 
Q in an interval J C [0, 1] can be expressed as 



#(Jn^ g ) 



QV\ 

2C(2) 



+ 0(QlnQ). 



Recall that if - < ^ are two consecutive elements in To, then 

a q — aq' 



1 and q + q > Q. 

Conversely, if q, q' G {1, . . . ,Q} and q + q' > Q, then there are a G {1, . . . , q — 1} and 
a' G {1, . . . , q' — 1} such that | < |j are consecutive elements in Tq. Proofs of these well- 
known properties of Farey fractions can be found for instance in |26| . |23| . |30| . 

Throughout the paper we shall denote by Tq, and respectively by Tq, the set of pairs 

(f ' f 7 ) °f consecutive elements in Tq with q < q' , and respectively with q > q' . We also set 



{(a,6)eZ 2 ;gcd(<i,i>) = l}; 



and ^2 



E 



{a/q,a'/q')eT> 



a/q {a/q,a'/q')eT< a/q 

a/q£j a/qdJ 

Aq = {(x, y) G Zp r ; < x, y < Q, x + y > Q}; 
Ttm,n = [ m , m + 1] x [n, n + 1], m, n G R. 
For each region TZ in M? and each C 1 function / : 1Z — > C, we denote 



oo,n = sup ll^/lloo,^ = sup 

(x, y )en (x,y)en 



dx 



x,y) 



+ 



dy 



x,y) 



The notation f <ti g means the same thing as / = 0{g); that is, there exists an absolute 
constant c > such that |/| < eg for all values of the variable under consideration. When 
the constant depends on a parameter 5, this dependence will be indicated by writing / g. 
The notation / x g will mean that / <g and g <C / simultaneously. 

We shall be mainly interested in consecutive Farey fractions | < ^ in .Fq with the property 
that, say, | belongs to a prescribed interval J C [0,1]. The equality a'g — a</ = 1 yields 
a = q — q' , where x denotes the unique integer in {1, 2, . . . , q — 1} for which xx = 1 (mod (/). 



Thus | G J = [ti , *2] is equivalent to q' G J q 



(i) 



[(1 - t 2 )q, (1 - ti)g]. Moreover, ^ G J is 
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(2) 

equivalent to q € J^, := [t\q ,t2q], where this time q denotes the multiplicative inverse of q 
(mod q'). 

An important device employed in pjj, |2|, HQ to estimate sums over primitive lattice points 
is the Weil type [SSI estimate 



(2.1) \S(m,n;q)\ <C t(o) gcd(m, n, q) 2 

on complete Kloosterman sums 



„, \ ( mx + nx \ 
S(m,n;q)= ^ e ( ~ J 



xe[l,q] 
gcd(x,q)=l 

in the presence of an integer albeit not necessarily prime modulus q, proved in |27| (see also 
[T!^. The bound from (|2.1j) can be used (see [4, Lemma 1.7]) to prove the estimate 



(2.2) N q (l,J) = |j| \j\ + O f (gi+*) 



for the number Nqil^J) of pairs of integers (x,y) G X x J for which icy = 1 (mod q), 
whenever X and are intervals which contain at most q integers. 

We shall use the following slight improvement of Corollary 1 and Lemma 8 in 4 a . The 
proof follows literally the reasoning from Lemmas 2, 3 and 8 in [3]. 

Lemma 2.1. Let £l C [1,R] x [1,R] be a convex region and let f be a C 1 function on fL 
Then 

(i) ^ a ' & ) = 77^ / / f{x,y)dxdy + £ R ,n,f, 

£ R ,n,f<.\\f\\oo,nRlnR+ £ p/IU^lni?. 

(m,n)eZ 2 

T^-m.n CZ 

(ii) For any interval JC [0, 1] one /ias 

X] /( a > 6 ) = TT^y / / y) rfx d y + ^R,a,f,J, 
(a,6)ennz2 r {/ 

u>/iere 

^fl,n,/,J «5 ||/||oo,nm / i?l +5 + ||/|| 0Oin length(50)lni?+ £ P/||oo,7W lnJ2 

(m,n)eZ 2 



/or any 0" > 0, where b denotes 3, the multiplicative inverse of b (mod a), J a is either or 

Ja 2 \ and mt is an up 
functions y 1— ► f(x,y). 



(2) 

Ja , and mf is an upper bound for the number of intervals of monotonicity of each of the 



When writing 6 G J a we implicitly assume that gcd(a, b) = 1. 
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The proof of (ii) relies on (|2.2jl . We also note the following important corollary of (|2.2j) . 
which will be often employed in this paper and in the subsequent work from and 

Lemma 2.2. Assume that q > 1 is an integer, X and J are intervals which contain at most 
q integers, and f : I x J — > R is a C 1 function. Then for any integer T > 1 one has 



E /(a, 6 ) = // -ft*' + 



afe=l(mod g) 



where 



E q ,r,j,f, T «s ?V + 1/l|oc + Tq^WDfU + |J| lJl l Dfllc 



for all 5 > 0. ifere || • ||oo denotes the L°°-norm on I x 

Proof. If T > g, then the error is larger than the sum to estimate and there is nothing to 
prove. 

If T < q, we partition the intervals I and J respectively into T intervals I\, . . . ,I T and 
Ji,...,J T of equal size |2j| = ^ and \ Jj\ = The idea is to approximate f(x,y) by 
a constant whenever (x,y) G Zj x For, we choose for each pair of indices a point 
( x ijiUij) G Xj x for which 



(2-3) JJ f = Pi\\J j \f(xi j ,y ij ) 

For (x, y) G Xj x J7j the mean value theorem gives 



f(x,y) = fi.r, r! ,,j) + 0((|X i | + UiDp/IU) 



(2.4) 

This gives in turn 



n.n.hiir^+0[^,\\Df\ 



E /M) = E E 

a<=:l,b£j i,j=l (x,y)eXiXjj 

(2.5) a6=l(mod q) xj/=l(mod q) 



Since each interval Z and J^- contains at most g integers, estimate Q2.2JI applies to them 
and gives 

(2.6) N q {Ii t Jj) = ^ \Jj\ + 5 (q^ +S ). 
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As a result of (j2.fi |) and (|2.Hj) . the main term in (|2.5|) becomes 

T 

^ £ |T J ||^-|/(x^,^)+0 5 (T 2 ^+ 5 



<7 2 



J f +0 5 (T 2 q^ S \\f\U, 



while the error term in ()2.5j) will be 



□ 



3. The second moment of the first exit time for cross-like scatterers 
Throughout this section we keep < e < i fixed, and take 



Q = Qs 



the integer part of — 

e 



We also denote 



z 2 * 


= Z 2 \{(0,0)}, 




C e 


= {0} x [-e,e] U [-£,£] x {0}, 




v £ 


= {0} x [s,e], 






= inf{r > ; (t cos a;, r sin a;) G C £ + 


Z 2 *} 


t P 


= the slope of the line OP, 




x,y)\\ 


= \J x 1 + y 2 , x, y G M. 





Let 

£ e = C7 £ + {(g,a); a/q£T Q } 

denote the translates of C e at all integer points with slope in J-q. 

For each point A(q, a) with ~ G ^"q we construct a vertical segment iVS of length 2e and 
a horizontal segment VF-E of length 2e, both centered at A. 

Performing symmetries with respect to the integer vertical and horizontal lines, the problem 
translates into a covering version in M 2 . It is clear that one can discard the points (q',a') 
with gcd(g',a') = d > 1, which are already hidden by (^-, ^-). 

The trajectory will now originate at O = (0, 0) and end when it reaches one of the compo- 
nents (q, a) + C £ of <£ e , | G J-q, as seen in the next elementary but useful lemma. 

Lemma 3.1. Any ray of direction uj G [0, j] which originates at O inevitably intersects <t £ . 
Moreover, if 7 = | < 7' = ^ are two consecutive Farey fractions in Tq and tanw G [7,7'], 
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then the ray of direction uj intersects either (q,a) + C £ or (q',a') + C £ and does not intersect 
any other component of <£ £ . 4 

Proof. We shall utilize the inequalities q + q'>Q + l>^>Q> maxjg, q'}, getting 

a a' — e a + e a' 

IA = ~ <ts' = ; — <tN = < M.' = —. ■ 

q q> q q' 





FIGURE 3. The case q' < q and ts' < tw, respectively q' < q and t$> > t\y. 



Equivalently, the intervals I a / q = l 2 ^, \ £ Fq% cover [0, 1] and two such intervals I a / q and I a i/ q i 

overlap if and only if - and ^ are consecutive elements in Tq. 
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In the case q < q' , we set {W } = OW n NS and {%} = ON n N'S' (see Figure El and 
note that a < a'), inferring that 

arctan 7' 

l 2 E {u) duj = 2 area(AOAiV) + 2 area(AOAfo^') - 2 area(AAW W) 

arctan 7 



In the case q > q' one has a' < a. Moreover, 

f a a 1 — e\ , 

t a = 7 < nimftv^,^;) = mm , ; — < max{tw,ts't < 7 = 

\q-e q' J 

This shows that any ray of slope tanw G [7,7'] intersects either (g, a) + C e or (q',a') + C £ 
and no other component of £ £ (see Figures E] and |3J) . 

Besides, we estimate the average of the second moment of the length l e (uj) of the trajectory 
when tan to G [7,7']. 

When t s > < t w (i-e. a' + q > e + §), we set {S } = OS" n AW, {S } = OS" n AT5, and 
note (see the first picture in Figure OJ) that 

arctan 7' 

J = 2 area(AOA , S') + 2 area(AOAS ) - 2 area(AAS S ) 

arctan 7 



eg +g( g , -a) +0(e 2 



g' 

When i 5 / > t w , we set {W } = OVF n NS, {S' } = OS' n iVS, and get (see the second 
picture in Figure 

arctan 7' 

y duj = 2 area(AOA , S') + 2 area(AOAS ) - 2 area(AAWW ) 



(3.3) 



arctan 7 



= g- £ (g 2 -g' 2 ) +Q(g2)- 
g 

□ 

We consider the region 

= {(*> v) g ^ 2 ; 1 < s < y < Q, x + y >Q} 

and the function 

/(*, y) = * + £ fr 2 -» 2 > = + sx, (x, y ) g ng. 
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(0,Q) (1,Q) (x,Q) 



(Q.Q) 




(0,0) (x,0) (0,0) 

Figure 4. The region Uq 



Consider also I = [a, (3] C [0, ?], take ii = tana, £2 = tan/3, and let J = [ii,^] Q [Ojl]- 
For (x, y) G S7q one has x>Q — y>^ — 1 — y, which gives 1 — ey < e(x + 1) < 2ex. It is 
also seen that 1 — ey > 1 — |4 > 0. As a result we find that [|/[|oo,fiq < 3. Since e 2 ^!Fq < 1, 
formulas (jSU), <|3~% (f3~3l provide 



arctan 



(3.4) f l 2 £ {u J )du J = 2 J Y J J l 2 £ (co)duj + 0(l) = 2 J Y J f(q,q')+0(l). 

I a ll arctan 2 a/<2 

9 



To master the latest sum, we aim to apply Lemma 12.11 to ^q. With the notation from 
Section 2, relation QM.4JI yields 



(3.5) 



J ll(u)du = 2 E /(<*, &) + 0(1). 



(a,b)en Q 



We also see that for (x, y) G Qq one has 
l 1-2 ^ < hence 



\d£\ 



v(i-ey) 



t— 1 < e + 2S£ < 1 and ,^ 

i — x — x I ay I 



I 9/ I 



(3.6) 



E P/ik^<E E ^ «Ei = q. 



(a,ft)eZ 2 



x=l y=max{Q— x,x} 



x=l 



Now we can apply Lemma 12.11 (ii) to the sum from Q3.5|) . and employ (|3.6j) and mj < 2, 
to infer that 



(3.7) 



/ I'M 



duj 



2{t 2 - h 
C(2) 



/(x,y)dx(iy + 5 (Qt+ <5 ). 



12 BOCA, GOLOGAN AND ZAHARESCU 

When a = and j3 = j, Lemma l2~Tl (i) improves upon the error in (|3.7|) to 

tt/4 

(3.8) J l 2 £ (u)du = -^-JJ f(x,y)dxdy + 0(Q\nQ). 

Q Q 
In summary, ()3.7|) . (|3.8j) and the equality 

J J f(x, y) dx dy = 1 2 & 

lead to 

tt/4 

Theorem 3.2. (i) e 2 J l*{u) duj = 1 + 2 n2 + q( £ | mg |) as e ^o+. 

o 

(ii) For any < a < (3 < f and <5 > ; one /ias 

9 /",9, x , (1 + 2 In 2)(tan /3 — tan a) „ , i x . , 
e 2 1 2 £ {uj) dw = H - + 5 (e2- s ) as e ^ H 

a 

Part (i) of this result was already proved in [221 • 



4. The r th moment of the first exit time for cross-like scatterers 

In this section we estimate the average of the first exit time for crossdike scatterers, 
thus proving Theorem 11.21 The first step towards estimating the integral Jjl r £ (<ju) dco = 
Jj ll~ 2 (u>)l 2 (lo) du> consists in approximating l r £ ~ 2 by a step function. 

We take / = [a,f3] C [0, f], t\ = tana, ti = tan/3, J = [ii,£2] ^ [0,1]- For consecutive 
Farey fractions f < y from J n !Fq, where Q = [~] , we denote 

a a + e . a' — e a' 

u)\ = arctan — , u>2 = arctan , uj 2 = arctan — , lo^ = arctan — . 

q q q' q' 

The function l r £ ~ 2 will be approximated by the constants l £ ~ 2 (uii) = ||(q,a)|| r_2 on [u^a^] 
and by Z^" 2 ^) = \\(q J , a')|| r ~ 2 on [0^2,^3] when q < q' , and respectively by 11~ 2 (ll>i) on 
[wi,^] and by l r £ ~ 2 (uj^) on [cj 2 ,W3] when q > q' . To be precise, we set 

UJ3 



A r ,j,e= J Eii(^ a )ir 2 / ifH<^+ J En (<?>') ir 2 / ^ 2 m^, 
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Next we estimate the quantities 



r,J,e 



and respectively 



E 



(2) 
,J,e 



U>3 

J r E (cj)du-A r>Jje 
^ J f r e (u)du-B rtJ)ts 

.111 J 



< E 



(i) 

r,[0,l],e 



< E 



(2) 

,[0,l],e 



o-'/q' ui 

An inspection of the case q < q 1 in the proof of Lemma 13. II leads to 



sup |/r 2 H < \\(q,a + e)\\ r -' - \\(q - e, a)||^ « r eQ 



r-2 



I r-2 



y-3 



CJG[ii!l,(i!2] 



and to 



sup IZr'H-CM \<\\W,a! 



I „i\\\r-2 



^£[1^2,1^3] 



, (q + e)g' 
9 > 



r-2 



, (a + e)q' 



Q r - 3 « 



r-3 



Therefore 



(4.1) 



W3 1^2 ^3 



< r q z (uj 2 - ui)eQ T 3 + q' 2 (uj 3 - uj 2 ) 



Q 



r-3 



But w 2 - wi < ^ - - = £ and w 3 - w 2 < ^7 - — = M 2 - < J r, so the right-hand side 

1 — g q q z — g' q qq' qq' ' 



dHJ) is 

As a result we infer that 
(4.2) 



a/ _ a+e 1— eq' ^ _J_ 

q ~ qq' 

r-2 r>r-2 



« r qe 2 Q r ~ z + 9L- « £ 



<?o,i] j£ = Or^- 2 E^=O r( Q^lnQ). 
In the case q' < q we get (in both subcases £5/ < i^y and is" > 



sup |c 2 M-CV)|< 

ue[uijj' 2 ] 



q, 



(a' - e)q 



r-2 



\\(Q,a-s)\\ 



r-2 



<Cr 



and 



sup - r 2 (^3)| < IK^.oOir^ - WW, a' - e)\r 2 «r eQ r 



ir-3 
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Employing also J 2 - w a = ^ - | = < j- q 
the case q' < q the estimate 



t 2 < 7^7 and w 3 - wi = < - *M = 4, we get in 



a' a' 



^3 



W3 



£M da; - / I'M du ~ ITM) / I'M du; 



qt-1 Qr-3 QT-2 



qq' 2 q' 



Hence 



r-2 



,/2 



■ (9,9')eA Q 



Since the contribution of a sine le term £f £(w)du; is « max{t? ; (? ' }r < Q^ 1 , we infer from 
(IP) and (gini) that 

(4.4) 



J l r £ {u)duj = S r ^ £ + O r {e 1 - r ). 



Next, we adjust the second integral in the expression of A r j £ , writing 



/2 , 12 12 
a + q = q 



7 + 1 = * 



/2 



'.: ii 



(a 2 + g 2 ) 



1 + 



0V + * 2 )(i + o 



This gives in turn 

(4-5) ii(9V)ir a 

In a similar way 
(4-6) \\(q,a)\r 2 - 



r-2 



-) \\(q,a)\r 2 (l + O r (- 



2 / 
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For further use, it is also worth to note 
{ 4 7) %— = — Q ' = q — 1 + O I — 

' } i + r 7 r (i + (f) 2 )(i+o(i)) i + (f) 2 ^ + \q 

Making use of (see the proof of Lemma 13. 1|) 

J l 2 £ (uj)duj = 2area(AOAiV) + 0{e 2 ) = sq + 0(e 2 ), 

UJ- A 

j l 2 (oo) dw = 2area(AOA^') = ' ^ ' 



q 



and 1)4.5(1 . we see that can be expressed as 

i\r-2 



(4.8) a/<? 



+ e ^||(g,a)|r 2 g + O r (l). 

If x denotes the inverse of the integer x (mod q) in [l,g], then a'q — aq' = 1 gives a = q — q'. 
Since < e, the error in the first sum in (|4.8|) is <C r Q 2 Q r ~ 1 eQ~ 1 = Q r-1 , and so ^4 r ,j,£ 

is equal up to an error term of order O r (Q r ~ 1 ) to 

Ir—X c-Jr 



Eii<«.«)ir , («»+£r-f&) 



a/q 

( 4 - 9 ) « „ . .. r-,/ T r-1 



9—1 max{5,Q— q}<a;<Q 



(1) 
J 

When q > q' , we employ (see the proof of Lemma 13. 1(1 



with J q as defined in Section 2 



w 2 

j l 2 £ {u) duj = 2 area(AOA^) + 0{e 2 ) = (1 - ^ Q + 0(e 2 ), 
J 1 2 (uj) duj = 2area(AO^'5') = eq' , 
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together with Q r ~ 2 £ 2 #J r Q < Q r ~ 2 , relation l|4.fij) . and the fact that a'q — aq' = 1 implies 
a! = q (mod q'), to infer that -Br,j,e is expressible as 

E J ll^«)ll r '- 2ii -^ + £ J (~Yli(9. *)irV + OrCCT" 1 ) 

a'/g' q a'/qi 

= IK*V)ir 2 (V + £ - + OrW- 1 ) 

a'/g' \ « « / 

= E E (^ 2 + ^^(^ + S-^ : r)+o r (Q^), 

g'=l max{ g ',Q- g '}<a;<Q V 7 

8 

where x denotes the inverse of an integer x (mod q') in [1,</]. Changing notation, -Br,j,£ can 
be rewritten as 

(4-10) E E (- 2 + 1 2 )^(^+^-^l)+Or(Q r ~ 1 ). 

9=1 max{q,Q-q}<x<Q \ V / 

*eJ< a > 

By and (jUTDJ), we infer that 

(4.11) y £(a/) ^ = T r ,J,e + O r {Q r - 1 ), 

l 

where T rj j(e) = 5i(Q) + S 2 (Q), with 

Q 

(4.12) S fe (Q) = ^ £ fk(x,x,q), k = 1,2, 

9=1 max{Q-g,g}<x<Q 

and 

f 2 (x,y,z) = (y 2 + z 2 )^ (ez + 

fi(x,y,z) = f 2 {x,z- y,z). 

For each q £ [1,Q], the functions /&(•, •,(?), defined on [1,(5] x [l,g], manifestly satisfy the 
estimates 

(4.13) || /fc ( Vjg )|| M<r 
and 

(4.14) p/^., . jg )|| oo « r .^l!<^_l. 

q q 2 

Thus we may consider in Lemma 12.21 for each q 6 [1,Q] the function /fc(',-,?), the intervals 
T = (max{Q — and J = J^p with \T\ < q and | J7"| = q\I\, and take T = [Q a ], to 



z' r 1 ex r 



z 7 - 1 z r ~ x 
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infer that the inner sum in (|4.12|) can be expressed as 

Q 

dx I dy fk(x,y,q) 



<p(q) 



I 2 

max{Q— q,q} j(*0 

Summing up over q £ [1, Q], we arrive at 

(4.15) s k (Q) = — g(q) + OsAQ r ^ +2a+5 + \i\Q r ~ a ), 

U q 



where 



Q (l-ti)z 
g(z) = ^- j dx j dy fi(x,y,z) 

max{Q— z,z} (1— £2)2 
i 2 2 



max{Q— z,z} tii 



The formulas 



/ Q bz \ Q bz 

~dz\~z dx d yh{x,y,z)\ =-— j dx j dyh(x,y,z) 



3. 

z az ' z az 



Q bz bz 

+ ~zj dX I dV lh^ X,y,Z " > ~~zJ h ( z,y,z ^ dy 

z az az 



H — / h(x,bz,z)dx / h(x,az,z)dx, 



(Q bz \ Q bz 

- J dx J dyh(x,y,z) \ = --^ J dx J dyh(x,y,z) 
Q~z az ' Q-z az 

bz bz 



1 f f dh 1 f 

+ - / dx dy — (x,y,z) + - h(Q-z,y,z)dy 



H — [ h(x,bz,z)dx / h(x,az,z)dx, 

z J z J 

Q-z Q-z 
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and the estimates (|4.13|) and 1)4.14)1 show that <C r — z — • As a result we get \g'(z)\ dz <C r 

Q T ~ 1 \a.Q. It is also clear that ||<7||oo Q r l , so we are in the position of being able to 
apply Lemma 2.3 in [3j to g, collecting 

Q Q / Q \ 

E *Y 5(9) = ^2)j 9( z ) dz + ol (ll^lloo + J \g'(z)\ dz^j In Q j 

9=1 1 \ 1 J 

Q 

= ^-jJg(z)dz + O r (Q r - 1 ln 2 Q). 
1 

Comparing the previous relation with (|4.15|) . we infer that both S\(Q) and ^(Q) can now 
be expressed as 

Q 

(4.16) -L J g(z) dz + OsAQ r ^ +2a+S + \I\Q r ~ a )- 

1 



Taking into account 1)4. and ()4.12jl . we gather 

Q 

(4.17) J r e (u) dw = JL J g ( z ) dz + OsAQ r ^ +2a+5 + \i\Q r ~ a )- 

I 1 

Integrating with respect to y in the formula that gives g and changing then z into Qz, we 
may express the main term in 1)4.17)1 as 

hz 1 dx ( £ + X —-^V-i 



(4.18) 



C(2)i J V * r 

1 max{Q— z,z} 

" ,J - f dz f dx ( ez r ~ x + 



C(2) J J \ Qz z 

l/Q max{l-z,z} 



where 



t 2 

K rI = 2 / (1 + t 2 )^dt = 2 



dx 



COS' X 



1 



Up to an error term of order O r (Q 2 ), the double integral in the right-hand side of 1)4.18)) is 
given by 

1/2 1/2 1/2 

e z(z r 1 + (l-z) r - 1 )dz + — / — ^ ^—dz — — / y '-—dz. 

J rQJ z(l-z) r + lj z(l - z) 

00 
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Since Q r+1 = -r+r + O r (pr), we now infer from IJ4.17J) and (j4.18j) the equality 

(4.19) J ^ = %J^- + 5>r (s- r+ ^- S + \I\e-^), 

i i 

with c r as in Theorem ll.il The proof of Theorem ll.2l is now complete. 



5. The moments of the number of reflections 

We take as before Q = [~] and keep up with the notation from the beginning of Section 4. 
An inspection of the proof of Lemma 13.11 shows that if | < ^ are consecutive in Tq , then 

if tana; G [- 



R £ (co) 



qi q-e' 
a a+e ' 



q + a 

q + a + 1 if tan w G [ 
q' + a' iftanw€[2±£,|] 
'g + a if tanw G [§, ^) 
y + o' if tanwe [^,^], 



+ a 



iftanwelf,^) 



+ a+1 if tanw G [— , ^ 

1 Lq— e' <r 



if g < g'; 
if g > and £5/ < tw; 
if g > q' and tg/ > tw- 



[q'+a' iftanoJG [^,§r], 



A first immediate remark is that we may replace q + a-\-\ by g + a in the above formulas, 
since the contribution of the corresponding arcs is small, as we see from | arctan x — arctan y\ < 
\x — y\, and from 



E 

a/q&FQ 



a + e 



Ee(q — a — e) 



and 



E 

a/q€T Q 



a — e 



q' q — e 
As a result we may write 



a/q^Q 

l-e(q + a') + e 2 



> - < 1 

a/q&^Q 



E 



i/q&pQ 



q'(q-e) 



« E A 



1. 



a/q£FQ 



where we set 



n(l) 



V~V<7 + a) r f arctan — : arctan — 

a/q 



r-l\ 



J v^/ / r\r ( a ' a + e\ 
+ > (q + all arctan — — arctan 

V V 1 1 J 

a/q 
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and 



T^ 2 } = (q + a) r [ arctan ; arctan — ] 

r,J,e y q , q 

a/q 

EJ . ( a' a' — e 
(q + a) arctan — - — arctan 
V i i 



a/q 

Employing 



h 



arctanfx + h) — arctan x = z + 0(h 2 ) = — ^ + 0(h 2 

1 + x z 1 + (x + h) z 

we now arrive at 

(5.1) J R r £ {uo) du = S r ^ e + T r ,j, £ + OriQ 1 - 1 ), 



where 
(5.2) 

and 
(5.3) 



l-eq' 



1—sq 



a/q 

To further simplify the expressions in Q5.2JI and (|5.3j) we employ 



a' a 1 a „ ( 1 \ / a\ ( / 1 



and 



to infer that 



9/ V K Q J J 



a/q&r Q a/q&FQ 



a/q { Q } ^ q J 



r-l> 



and also that 
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The main term in (|5.4|) can now be conveniently expressed as 



(5.5) 



.4 



r.J.e 



5=1 max{Q— q,q}<x<Q 



xe J, 



a) 



where x is the multiplicative inverse of a; in {1, — 1}, and this time we set 



f(x,y) = f{x,y,q) 



(1 + gr 

g2 + ( ™ 



£q r-l + L^ x r-l 



q(2q — y) r ( 1 — ex / x \ <- J 



g 2 + (g - y) 2 

Since < 1 — ex < 1, it is easy to check that 

(2g - yf ^ ,-2 



9/ 



dx 
Of 



oo,Xx J, 



g ( 2g - y r j_ 3 

(i) g + (g — y) g 



<9y 



oo,Xx J, 



(2) 



9 / (2g-y)' 



<r g 



r-3 



dy \q 2 + (q- y) 2 , 

Taking these estimates into account and applying Lemma \'2.2\ for each q G [1,(5] to / = 

/(•, -,q),I = (max{Q — (/, (/}, Q], i7 = Jq 1 ^ and T = [Q a ], we infer that the inner sum in (|5.5() 
can be expressed as 



e + 



1 — ex /x\ r ~ 1 



m&x{Q-q,q} 



q V <7 



dx 



r(U 



g(2g - y) r 
g 2 + (q - yf 



dy 



+ O r>5 (Q 2a q^ +S q r - 2 + Q a ql +5 q r - 3 + |I|<rV? r ~ 3 ) 

g r ~V(g) + O r , 5 (Q-f +M + |T]Q— ) 5 



where 



and 



C, 



tT-1 



(2g - y) r 



dy 



q 2 + (g - y) 2 g r 1 7 q 2 +y 2 



t2q 

j_ Uq + y y 
r-i y - 



dy 



(l-ta)g 



*i<2 



*2 

y ill ^ df = J(l + tanx) r dx 



ti 



Q 



9r(q) 



e + 



1 — ex /x\ r ~ 1 



m&x{Q—q,q} 



q v g 



fix. 
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We now arrive at 

(5.6) V f{x, x,q) = Q± h r (q) + O r ^Q r -'^ 2a+S + W^ 1 ) 

^-^ q 

ma,x{Q—q,q}<x<Q 



h r (q) = I I eq 1 " 1 + ^— ^ x r ~ l ) dx. 



where 



max{Q— q,q} 

Since ||/i||oo <C Q r_1 and f® \h'(q)\ dq < Q r ~\ Lemma 2.3 in together with (J53J) and 
(|5.6|) show that 

9=1 q 



^ J h r (q) dq + O r , 5 (Q r ~5+ 2a + 5 + \I\Q r ~ a ) 



C(2) 
1 

Making use of = 1 + 0{e) we arrive by a straightforward computation to 
<2 l 



J h r {q)dq = Q r j (^x r 1 (l — max{l — x, a;}) + 



1 — max{l — x, x} 1 



rx 

_ l-max{l-,,,r +1 ^ x + 0r(Q ,- 1) 
The integral above is seen to coincide with 

1/2 

r-1 ,1 xr-U 1 " (1 " X) r 1 " (1 " X) r+1 \ , TT 2 C r 

x[x r 1 + 1 - x) r x ) + £ f- - v . ./ — - dx - 



rx(l — x) (r + l)x(l — x) J 12 



o 

hence 

o 



j K(q) dq = TLZ9L + CW" 1 ) = ILSdEL + O r (e~^) 
l 



and as a result 



S r ,J,e = 4 • "^Sr- A 1 + tan2 r + O^^e-^^ 2 "" 5 + |/|^+") 
= J (1 + tan x) r + O r . i5 (£- r+ ^ 2a - <5 + |/|e- r+a ). 
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By reversing the roles of q and q' it is seen in a similar way that 

„ c r e~ 



2H 



-r,J,e 



(1 + tan x) r dx + O r {e- r+ ^- 2a - 5 + \I\e~ 



This concludes the estimates of S r j e an d T Ty j )£ . Theorem 11.51 now follows from ()5.1|) , 

6. The case of circular scatterers 

Note first that the statements of Theorems 11.21 and 1 1 . 51 hold true if we replace the scatterers 
C £ + Z 2 * by V e + I?*. 

In this section we consider the circular scatterers D £ + 7?*, where 

D £ = {(x,y)£R 2 ;x 2 + y 2 = e 2 }. 

For each integer lattice point (q, a), let (q, a±e±) denote the intersections of the line x = q 
with the tangents from O to the circle 



D 



e,q,a 



(q, a)+D e = {{x, y) G R 2 ; (x - q) 2 + (y - a) 2 = e 2 }, 



where e± = £±{q,a) are computed from the equality 

£ = 



which gives in turn 



or 



e±q 



l + (^±) 2 V / q 2 + (a±£±) 2 ' 

e 2 ± q 2 = e 2 q 2 + e 2 (a±£±) 2 , 
(q 2 - e 2 )4 T 2ae 2 e± - e 2 (q 2 + a 2 ) = 0. 



The latter provides 




Figure 5. A circular scatterer 



ae 



(6.1) e± = ± ^p 9 + 9 ^ ~ y/<f+a*q* -a*e*. 

qZ _ £ z gZ _ £ z 
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Employing also 



vY + a 2 q 2 - vY + a 2 q 2 - a 2 e 2 < e 2 , 

2, 

q z — e z q z 



and 



we arrive at 



ae 



q 2 — e 2 



(6.2) 



q J cos arctan - 



Proof of Theorem ll.il We wish to compare L f £ (u>) duj with L f £ (tu) duo where f £ (co), the 
smallest r > for which 

(r cos lo, t sin to) G |^J {q} X [a — £-(<?, a), a + £+(q, a)], 

denotes the first exit time in the case where the scatterers are the vertical segments {q} x 
[a — e-(q, a), a + e + (q, a)]. From Figure it is apparent that 

sup \f E (u) - f e (uj)\ < 2e, 

LU 

and so, since sup^ t £ (uj) < sup w 1 £ (uj) < we get 

sup|fJ(o;) - tJ(cj)| < r e 

which gives 



V2 



r-l 



(6.3) 



f r E {u) dcu = / f r £ {u) dw + O r (e 2 - r ). 



To estimate fjf £ (uj) du, we divide the interval / into N = [e e ] intervals of equal size 

in 

N 



Ij = [ujj,ujj + i] with \Ij\ = U x e 9 for some < 9 < |. Then one has for all j that 

COSCJj COSWj + l 



-i 



thus the integers = [ 

(6.4) 
and 

(6.5) 



+ 1 and Qj = [ ^a/V ] satisfy 
< Qt - QJ « s 6 - 1 



1 e _ £ 3/2 £ + £ 3/2 1 

— r < < < — 

cosuj cosw i+ i Q 
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Furthermore, it follows from (|6.2j) that there exists Eo = Eo(6) > such that for all e < £q, 
all j, and all | S [tancjj, tancjj+i], one has 

e-e 3 / 2 . . e + e 3 / 2 
< e±(q,a) < . 

COS UJj COS UJj + i 

This implies in conjunction with (|6.5|) . for all - G [tancjj,tanu>j + i], the inequalities 
(6.6) -L< s±(q , a) <-L. 

Since Q± = + O^" 1 ), one has 



(6.7) {Qt) r = ^f^ + O r [E-^) 

The first exit time increases when all the sizes of scatterers decrease. Thus we infer from 
(|6,6|) the inequalities 

(6.8) / V^{u))dw< [ f r e {uj)duj< [ l\ (Lu)duj. 

J QT J J Of 

I 3 T T ■ 3 



But by Theorem 11.21 and by Q6.7|) we may write 

(6.9) I V\ (u) dw = C r (Qfy f + O r , S (E- 



r+i-2a-<5 _|_ £ -r+6+a\ 



T 3 
1 1 



with the better error term e 2 s for r = 2. Also using [ T dx r <C r l-ZV I <C e e we infer that 

j COS 3^ i J i 

the first integral in ()6.9|) is expressible as 

u; 3 + l 

( 6 .10) CrCOST ^ [ J^ + 0rS (E~ r+2e + E- r+1 2-^- S +E- r+e+a ), 

£ r J COS r X 

with the better error term E~ 2+e + e"2 _<5 for r = 2. 

Summing up over j we infer from (|6.1U|) and (|6.3|) that 



TV , 
ax 



(6.11) 



y (w) = ^ ^ cos r y 



COS' X 



+ O ri5 (e- r+ 5- 2Q - e - 5 + e - r+a + E-' r+e ) 



with the better error term e~2~ e ~ s + e~ 2+e for r = 2. 

Finally, we apply the mean value theorem and chose some £j £ [u)j,Uj+i] to evaluate the 
sum 

N 



^2 cos r y 

i =1 



dx 
cos r x 
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as 



N , x r N 

ElWj+l — OJj) COS UJj v-^/ \(- ^ I \\ 

3=1 J 3=1 

iV 

= ^(^ +1 -^)(l + O r (e e 
= |/|+O r (e e ). 

This implies Theorem II .11 in conjunction with by taking 9 = a = | for r / 2 and 6 = \ 

for r = 2. □ 

Proof of Theorem \1.4\ We proceed along the same line to estimate the moments of i?. Here 

we denote by R(u) the number of reflections in the side cushions in the case of vertical 
scatterers (of variable size) {q} x [a — e-(q,a),a + e + (q,a)], (q,a) € Z 2 *. It is seen as in 

the proof of Theorem 11.11 that L R t £ {oj) du> differs from L R r £ (uj) du> by an error term of order 
O r (s 2 ~ r ). One can also show that 

(6.12) f R r _^(u)dw < [ R r E (u)du) < [ R r 1 (u)du>. 

J Q- J J Of 

J 3 J T 3 

2 J 1 ] 1 J 



Applying now Theorem EH to the vertical scatterers V\/q± on the intervals Ij = [uj,Uj+{\ of 
equal size = jj- xes with 6 = a = |, and also using (|6.7jl . we find that 

Ir\ (uj)duj = CrC °^ u i [ {l + tmixf dx + O r ^e- r+ \- 5 ), 
J of £ T J 

T 3 1 

and thus 

N u *+* 

(6.13) / R r £ {uj)duj = ^J^coe r Wj- / (l+tanx) r dx + O r)S (e- r+ *- 5 ). 
I 3=1 ujj 

By the mean value theorem we find £j , rjj E Ij such that 

N N 

(6.14) y~] cos r Uj / (1 + tanx) r dx = \^(uj+i — u)j) cos r + tangj) r , 
and respectively 

(sin x + cos x) r dx = / (sinx + cos a;) r tix 

/ 

iV AT 
= — Wj)(sinr/j + coS77j) r = ^^(uj+i — Uj) cos r rjj(l + tanr/j)' r . 

3=1 3=1 
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From 

cos r uij = cos r r/j + O r (ojj + i — uij) = cos r rjj + O r (e») 

and 

(1 + tan^-) r = (1 + tanrij) 1 " + O r (\ tanfj - tan^|) = (1 + twqjY' + O r (es) 
we infer that the sum in (|6.14|) is equal to 

N 

y^X^j+l - Uj)(^cos r rij + O r (e^)^j ((1 + tan^) r + O r (ei)) 
1 

This can be combined with 13|) to collect 

j R r e {oj) duo = J (sinx + cosxY dx + O r: s(e~ r+ ^~ S ), 
I I 

which concludes the proof of Theorem 11,41 □ 
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